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Within the macroscopic Landau-de Gennes approach, we examine the Gaussian normal mode 
fluctuations of semi-infinite nematic liquid crystals in a regime of critical wetting. It is argued that 
surface free-energy potentials which strongly suppress the long-range nematic order favor the ap- 
pearance of bound biaxial nematic-director fluctuation modes, located in the domain occupied by 
■ the thermodynamic phase wetting the wall. Instead, substrates enhancing the orientational order 

Q promote the existence of uniaxial nematic-director local excitations. Close to the phase coexistence 
temperature both types of local exciations are strongly softened as compared to their bulk coun- 
l^-} ' terparts and acquire characteristic cusplike low-energy spectra. These spectrum peculiarities are 

directly connected to the critical behavior of the mean-field interface position and can provide a 
valuable insight on the nature of surface interactions and critical wetting phenomena in nematic 
liquid crystals. Possible changes in the local director mode properties resulting from the critical 
interface position fluctuations and order electricity effects are also discussed. 



PACS: 61.30.Cz, 64.70. Md 



I. INTRODUCTION 
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Liquid crystals in confined geometries (such as plates, droplets, porous glasses, and polymer networks) constitute an 
interesting class of thermodynamic systems both from experimental and theoretical viewpoints In particular, the 

generic critical wetting phenomena |^||, studied in a large number of scalar order-parameter systems, are expected to 
acquire some new interesting feartures near the bulk nematic-isotropic coexistence temperature Tjvj arising from the 
complicated nematic order parameter (a symmetric traceless second-rank tensor) . In confined nematic liquid crystals 
the critical wetting transition usually appears as a continuous interface derealization transition at bulk coexistence, 
when the nematic-isotropic interface critically intervenes into the bulk as Tjvj is approached. Typical features of a 
critical interface derealization transition are the diverging interface position and the existence of a soft fluctuation 
mode directly connected to the interface position fluctuations (for a review see, e.g., Ref. §). At a mean-field level 
Q\ , and in a uniform-director approximation, the macroscopic Landau-de Gennes approach || reproduces the well-known 
■ critical wetting results originally established for scalar order-parameter systems, so that many known results can 
^2 ! be directly applied to nematic liquid crystals (for a review on wetting transitions in liquid crystals see Ref. § and 
p ■ references therein). 

Concerning the fluctuation dynamics, the situation in systems with an orientational rong-range order might be 
substantionally changed. Generally, the nematic order-parameter fluctuations are characterized by five independent 
fluctuation fields, one of them being related to the scalar order-parameter fluctuations, whereas the other four fields 
, describe the transverse fluctuations of the nematic director. Therefore, an interesting issue is to examine nematic 
director fluctuations in a regime of critical wetting, as one can expect that some peculiarities of normal-mode director 
. \ fluctuations will reflect the criticalities. Published experimental works show quite different elementary excitation 
spectra, ranging from very small to very large relaxation rates as compared to bulk excitation spectra [ p~0|JTl| ] . Up 
to now, relatively small number of theoretical studies treating the fluctuation dymamics in confined liquid crystals 
near the bulk coexistence have been published. Recently, Ziherl and Zumer jl2| have published a work dealing with 
order-parameter pretransitional fluctuations above Tjvi under the condition of a strong homeotropic anchoring of 
the nematic director. Apart from the typical soft fluctuation mode, reflecting the interface position fluctuations, 
their numerical analysis also indicated slow uniaxial director modes located close to the substrate. In the discussed 
work the critical wetting regime is ensured by a strong scalar order-parameter anchoring at the surface, describing 
substrates enhancing the nematic order. Another extreme type of substrate anchoring yet allowing for a critical wetting 
regime, called here weak- anchoring limit, corresponds to substrates strongly suppressing the nematic arrangement. 
For example, substrates treated by the S^O-evaporation technique may produce such disorder effects. In this case 
(and under some additional conditions discussed below) , it is the isotropic phase which wets the wall, so that one may 
expect different behavior of the normal fluctuation modes. 
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A goal of the present work is to compare the fluctuation mode dynamics based on such extremely different substrate 
potentials. The emphasis is especially on substrates strongly suppressing the nematic order. A characteristic feature 
of the weak anchoring is the possible critical behavior of the surface scalar order parameter by itself. This may also 
be of importance in the critical wetting experiments, as the finite size effects frequently mask the weak logarithmic 
divergence in the interface position (characteristic for short-ranged substrate potentials) As argued below, in the 
weak-anchoring limit, when the isotropic phase wets completely the substrate, in addition to the scalar order-parameter 
soft excitation, gapped biaxial director modes located in the presurface quasi-isotropic layer can exist. Basically, the 
low-energy spectrum of these local biaxial excitations is determined by the mean-field interface position. Since the 
latter is logarithmically divergent at T = T/vj, the excitation spectrum shows a characteristic cusplike behavior in 
that limit, as found also numerically in the strong-achoring case |T^ |. 

The cited results are based on the Gaussian approximation and suggest constant nematic-director mean-field con- 
figurations, when the director fluctuations decouple from the "dangerous" scalar order-parameter modes. In principle, 
a series of effects can invalidate the above two approximations. At a first place, higher-order fluctuation effects may 
be of importance since D = 3 is the upper critical spatial dimension for a critical wetting transition. In particular, 
it is well-known that the effective nematic-isotropic interfacial width diverges in D = 3 jl^] . Clearly, this latter 
divergence can modify the local excitation spectra. However, close enough to 7jvj the simple Gaussian picture of local 
director excitations is still valid, since the effective interfacial width is much smaller than the domain thickness (see 
below). As to the inhomogeneous director configurations, they can result from many sources, the main effect being 
the appearence of intermode couplings. As an example, we will discuss below order-electricity effects, as the problem 
is characterized by strong scalar order-parameter variations in the interfacial region. 

The rest of the paper is organized as follows. In Section II we define the surface Landau-de Gennes free-energy 
functional and the free-energy substrate potentials describing the strong- and weak-anchoring limits. For uniform- 
director configurations, we present the mean-field scalar order-parameter profiles which are later used to fix the relevant 
areas of the phase diagrams. In Section III we study the Gaussian order-parameter fluctuations. The emphasis is 
on the properties of the local nematic-director excitations in the two cases of substrate anchorings and in complete 
wetting regimes. The role of interface position fluctuations and order electricity effects are also analyzed here. The 
last Section contains discussions and a summary of the results. 



II. LANDAU-DE GENNES MODEL 



We adopt the macroscopic Landau-de Gennes description (one elastic-constant approximation) based on the fol- 
lowing free-energy density || : 

/lg(Q) = -^(d a Qf3y) 2 + -QaffQafi ~ T^QapQ /3yQya + J {Qaf}Qaf}) 2 ■ (1) 

Here Q = Q a p{^) denots the nematic order-parameter field (symmetric traceless second-rank tensor), L is an elastic 
constant, and a = a(T — T*). The material constants a, b, and c are supposed to be positive and temperature- 
independent. T* is the supercooling temperature for the bulk isotropic phase. 

It is convenient to introdice the following reduced quantities: (i) the reduced tensor order parameter S = Q/Q c , 
where Q c = 26/3 v^c is the value of the scalar order parameter Q at the bulk nematic-isotropic phase-transition 
temperature T NI = T* + b 2 /27ac, (ii) the reduced temperature r = (T - T NI )/(T NI - T*). The liquid crystal is 
assumed to occupy the semi-infinite space z > 0, z being the coordinate normal to the surface. Using the above 
reduced quantities, the surface free-energy functional can be recasted to the following form: 

F S [S}=LQ 2 C J dxdy J™ dz|i(a Q 5 /37 ) 2 + ^[/(S)-/(S 6 )]+5(z)/ (S)|, (2) 

where the uniform free-energy potential /(S) reads 

/(S) = i±Itr S 2 - V6tr S 3 + l - [tr S 2 ] 2 . (3) 

Here S& = Sf,(r) denotes the order parameter deeply in the bulk (z — > oo). The correlation length of the isotropic 
phase at coexistence is denoted by £o — \/ L/ oi(Tni — T*). The short-ranged substrate potential /o(S) (defined at 
z = 0) comes from a variation of the intermolecular forces near the surface. Up to second order in S, /o(S) is frequently 
modeled by the expression jl^] 
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/o(S) 



-tr(S - S 



2 

s J } 



(4) 



where S s is a symmetric traceless second-rank tensor characterizing the substrate, and a s is a phenomenological 
constant related to the ratio of the coupling constants at the surface to those in the bulk. 
The Euler-Lagrange equation related to Eq. (0) reads 



£o AS + (1 + r)S + 3\/6 



S 2 -i(trS 2 )l 

3 v ; 



2(trS 2 )S = 0, 



where 1 is the unit second-rank tensor. The boundary condition at z = 0, related to the last equation, is 

ds , a/ (s)" 



dz 



OS 



= 0. 



(5) 



(6) 



A. Uniform-director order-parameter profiles 

In a uniaxial uniform-director mean-field configuration, the solution of Eq. (|b|) has the following generic form: 

S ( z ) = -M(3no®no-l), (7) 
V6 

where no is the constant nematic director and u(z) is the scalar order-parameter profile. As discussed below, the 
uniform-director ansatz (0) is presumably not valid when Eq. (|2|) contains terms imposing elastic nematic-director 
distorsions (such as order-electricity or flexoelectric free-energy terms). Nevertheless, the uniform-director approxi- 
mation is a useful starting point even in the presence of such terms, as far as the system is in a regime of critical 
wetting and close to T/yj. 

For future references, we briefly study below the mean-field solutions of Eq. <M) as well as the surface phase diagrams 
(for the two extreme cases of surface free-energy potentials discussed above) in a vicinity of the bulk transition 
temperature T^j. In a uniform-director approximation, Eq. (Q) reduces to the following non- linear differential 
equation for the profile function u(z): 

d 2 u ldf(u)_ 



/(«) = - v? + \u\ (9) 



supplemented by the boundary condition 



du 

dz 



Here uq is the value of the scalar order parameter u(z) at the surface, uq = u(z) z =o- Deeply in the bulk the conditions 
read (du/dz) z=+oo = and u(z = +oo) = u^, where = u&(t) is the bulk scalar order parameter: ui, = in the 
isotropic phase and 



Ub(r) = - (3 + 8t) , t<0, (11) 

in the ordered nematic phase. 

It is instructive to rewrite Eq. (^), by use of the boundary conditions for u(z), in the form 

du 



^Y z = Sign(ub ~ «o)a/2/(u)-2/(u 6 ), (12) 

where sign(itb — uq) = ±1 for Ub > uq and Ub < Uo, respectively. Now, using Eq. (p^|), the boundary condition ( p^| ) 
can be transformed to an implicit equation for uq, 
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sign(u b - u ) — ^2f(u ) - 2f(u b ) = -—, (13) 
Co ouq 



and the free-energy functional (fq) reduces to the following function of u : 

fs(uo) 



F s [u ] sign(u 6 -u ) f /~oJT~\ o77 — V, / c \ d/A 

r ^, 2 = —, / du^2f{u)-2f{u b ) + f {u ). (14) 

LQc Co J„ 

In terms of the function J s (uq), Eq. ( |l3| ) is equivalent to the minimum condition df s (uo)/duo = 0, whereas the phase 
stability condition reads 

d 2 f s (u ) _ sign(u -M b ) df(u a ) | d 2 f Q (u ) ^ Q 
du 2 ~ S Qy /2f{u)-2f(u b ) du Q du 2 Q 

The last three equations completely describe the surface phase diagram of a semi-infinite nematic liquid crystal 
in uniform-director uniaxial states. In the context of surface phase transitions in nematic liquid crystals, using the 
substrate potential /q(i*o) = —hoUo, these equations were originally analyzed by Sheng [ fL6| . Surface free-energy 
potentials of the form 

/o(wo) = -hoUQ + yltQ, (16) 

constitute the basis of the theory of wetting fenonomena in simple liquids |5| and they have also been used in the 
context of liquid crystals by Sluckin and Poniewierski |17j . Since the emphasis in the present study is on fluctuation 
effects in a regime of critical wetting, we do not analyze the full phase diagram related to Eq. (||). Instead, we mainly 
consider uniaxial surface potentials which can be described by the tensor parameter 

S s = -^(3n s ®n s - 1) (17) 
V6 

in Eq. (0), where u s and n s are the scalar order parameter and the director orientation preffered at the surface. 
To be concrete, in what follows we assume that the unit vector n s is normal to the surface (homeotropic director 
anchoring). In the cases without elastic distorsion forces, such as those described by the free-energy potential (|l|) 
under appropriate boundary conditions, the director no is always parallel to n s , provided that u s > 0. In this case 
Eq. (^) is reduced to the following substrate potential: 

fo{u ) = y( w o ~ u s f, a s >0. (18) 

There are two interesting limits of Eq. (|l^) when the solution of Eq. ( |l3| ) is extremely simple: (i) The srong- 
anchoring limit, a s — > oo, which is appropriate for substrates imposing a strong scalar order-parameter anchoring at 
the surface. In this case the solution of Eq. ( [l3| ) just reads uq = u s . This is also the case discussed in Ref. [ ft2|| . 
(ii) The second interesting limit of Eq. ( |l8| ) corresponds to the choice u s = 0. This latter condition may fit to, say, 
substrates treated by SiO- evaporation technique, when the nematic order is strongly suppressed. In this case Eq. (|4j) 
reduces to /o(S) = (a s /2)trS 2 , so that the substrate does not impose any preferred direction on the nematic director. 
We call this weak- anchoring limit. It is important that the above simplified potentials generate the main critical 
wetting regimes found with the more general form (|l^), so that we restrict our interest to the above extremely 
different cases. 

For the bulk free-energy potential @ and for uniform-director mean-field configurations, the non-linear differential 
equation ( [l2| ) can be solved both above and below Tni and for arbitrary uq- The surface scalar order parameter uo 
should be independently obtained from Eqs. ( |l3|) and (|l^). The explicit expression for the order-parameter profile 
above T^i reads 

u(z)= 1 + T /, I v T>0 ' ( 19 ) 

1 + V r sinh + ai) 



cti = arcsinh 



1 fl + T 



Uq 



- 1 



where £/ = £o/\Al + r is the correlation length in the isotropic phase. 
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Below the bulk coexistence temperature Tjv/ the profile function reads 

3u b -2(l + r) 



u(z) = u b 



2u b - 1 - y/u h — 1 sinh [sign(u — u b )z/£ N + c(n} ' 



t < 0, 



(20) 



ajv = arcsinh 



V u b - 1 



2w b - 1 



3u fc -2(l + r) 

Uq - U b 



where £at = d 2 f (u b )/ du 2 = £,o/\/3u b — 2(1 + r) is the correlation length in the nematic phase. For r > the profile 
m(z) describes a paranematic layer, whereas the bulk is occupied by the isotropic phase. For large enough uq, u{z) 
has an inflection point at z = dj marking the center of the nematic-isotropic interface. The implicit equation for di 
is u(di) — ui, where the function u\ = Ui(r) = (3 — \/l — 8r)/4 corresponds to the local maximum of /(it), Eq. (|]). 
Explicitly, 



rf/ = £/ In 



u \ T + (1 - ui) + ^(1 + r) [t + (1 - ui) 2 



«i/ r+(l-« ) + V(r + r) [t + (1-u ) 2 



Uq> U\. 



(21) 



Similar formulae can be established below Tni- For uq < u b the surface layer is less ordered as compared to the 
bulk. If in addition uo < u\, the profile has an inflection point at z = d^ marking the interface position difined by 
the equation u(djv) = ui (see Fig. 1). The result reads 



dN = £/v In 



u b -u Q \ uiu b + (u b - l)(2u b + Ul) + yj u b (Au b - 3) [u 6 - 1 + {u b + ui- l) 2 



u b -uij Uo u b + (u b - l)(2u b + uo) + y/u b (4:U b -3)[u b -l + {u b + u - I) 2 ] 



(22) 



B. Surface phase diagrams near Tni 



The profile functions, Eqs. ([19|) and (p0[), may be used to describe the surface phase diagrams. In the context of 
interface derealization transitions [pl| , the surface phase diagram has already been analyzed in the framework of the 
substrate potential ( |l6|) @ . For our study it is enough to consider the extreme cases specified above and to mark the 
respective peculiarities connected to the form of the used surface free-energy potentials. 

(i) Strong- anchoring limit: 



Firstly, one must find the surface order parameter uo from Eq. (p_3[). In the limit a s — >• +oo the solution of Eq. (13) just 
reads uo = u s , so that the surface is characterized by the parameter u s . Let us consider the phase diagram for r > 0. 
If u s > 1 and r — > + , Eq. ( |2l| ) predicts a logarithmically diverging interface position given by di = £o ln(l/r) + 0(1). 
Exactly at u s = 1 one finds di = (Co/2) ln(l/r) + 0(l). The critical behavior is connected to the following non-analytic 
terms in the surface free energy (pf) considered as a function of r and A = u s — 1: 

fr g (r, A) = T - In A + 7grT7 + ^ v/A^ + i(A 2 + rfl\ (23) 
z r I 3 

In the region A < there are no singular terms at all, and the thickness of the presurface nematic layer di is of 
the order of the correlation length (partial wetting). On the other hand, one finds f^ m 9 oc rln(l/r) as A > 0. 
Thus, there is a line of critical points, defined by r = + and u s > 1, describing a critical wetting of the substrate 
by the nematic phase. The transition at r = + for u s > 1 is known as a complete wetting or as a critical interface 
derealization transition of type P + . (?]] The point (r, A) = (0 + ,0) is a special point on the phase diagram since the 
parameter A becomes a relevant field (in addition to r): the derivative d 2 f s /dA 2 oc A/V r + A 2 has a finite jump 
at A = 0. Respectively, the interface position di = £o m (l/|A|) + 0(1) is finite as A < 0, whereas it is delocalazed 
(infinite) when A > 0. As a matter of fact, the multicritical point (r, A) = (0 + , 0) corresponds to the tricritical point 
T + found for the substrate potential (|l^), but in the strong- anchoring limit (a s — > oo) the summetry breaking field 
ho is irrelevant. 

Similar critical properties can be established just below the bulk coexistence temperature Tjvj- Now the interface is 
located at z = djv, and the region < z < d^ is occupied by a suppressed (as compared to the bulk) nematic phase, 
provided that u s < u\. For positive values of the parameter u s , remains finite at r — 0~, and the thickness of the 
surface layer dN is of the order of £n- Complete wetting of the substrate by the isotropic phase is possible only for 
non-positive values of u s |llj. The point (r, u s ) = (0~,0) is a multicritical point analogous to those discussed above. 
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The interface position shows the same asymptotic behavior as di does: for u s < and r — > Q _ , it diverges as 
<iv = £o m (l/|r|) + 0(1), whereas at u s — one has djv = (Co/2) ln(l/|r|) +0(1). For fixed r = 0", djy is finite when 
?i s > 0, c?jv = Co hi(l/u s ) + 0(u s ), and it becomes infinite as u s < 0. 
(ii) Weak- anchoring limit: 

This case is more interesting in a sense that now the surface order parameter uo is not pinned by the surface and 
may change with the temperature, uq — uq{t). As obtained from the inplicit Eq. (|l3|), the function Uq(t) has the 
following asymptotic behavior. For large enough a s (a s > CcT 1 ) an( i r ~~ * — ) u o = ^o(t) vanishes as 

u (r) = -=^=|T|*+O(|r|) ) a^o > 1, r < 0, (24) 

whereas uo = all over the region r > 0. This additional surface criticality is specific for the weak-anchoring limit 
and takes place simultaneously with the critical interface derealization discussed above. The asymptotic behavior of 
dN as t — > 0~ (see Fig. 1) is analogous to those found in the case of srong anchoring: d^ = Co hi(l/uo) + 0(1), with 
Uo = uq(t) from Eq. (p4|). 

The case a s < C^ 1 ^ s ^ ess interesting in a sense that now the surface order parameter it is discontinuous at T^j: 
it jumps from the finite value uq = 1 — a s Co + 0(|r|) (for r — > - ) to zero (for r > 0). Respectively, there is no 
critical behavior for a s £o < 1, and d n is of the order of the correlation length in the nematic phase £/v- The point 
(t, a s £o) = (0~, 1) is a milticritical point (type T~ according to Ref. ( 0)) characterized by three relevant fields, i.e., 
r, a s £o, and ho- Exactly at a s £o = 1, the surface order parameter uq(t) acquires the asymptotic behavior 

i 

Mo(r) Hir +o(|r|) ' a >& =i > r<o - ( 25 ) 

Respectively, the interface position delocalizes as djv = Co hi(l/«o) + 0(1)), with no = uq(t) defined by Eq. fl25|), 

III. LOCAL FLUCTUATION MODES 

Within the Landau-de Gennes approach the normal fluctuation modes are obtained by solving the eigenvalue 
problem connected to the linearized Euler-Lagrange equation. The latter equation can be derived from Eq. (^|) by 
use of the order-parameter decomposition 

S(r) = S o (*) + 0(r). (26) 

Here So(z) is the mean-field solution (^), as defined by the scalar order parameter u(z), Eqs. ( |l9| ) and (|2p|). <fi (r) is 
the order-parameter fluctuation field (symmetric traceless second-rank tensor) . 
It is useful to adopt the following parametrization of <f>(r) [ po| : 

4 

0(r) = (27) 

i=0 

where the base tensors gi, i = 0, • • • , 4 read 

1 , 

go = -^(3n ® n - 1), 

V6 

gi = -^=(ei (g) n + n <g) ei), g 2 = -^=(e 2 <g) n + n <8> e 2 ), (28) 

g 3 = —^i e i ® ei - e 2 ® e 2 ), g 4 = ^|( e i ® e 2 + e 2 <8> e x ), 

and satisfy the orthogonality relation Tr(gi.gj) = The unit vectors ei_Le 2 are perpendicular to the nematic 
director n . 

The variable </>o( r ) describes longitudinal scalar order-parameter fluctuations, whereas the pair variables 
[0i(r), </> 2 (r)] and [03 (r), </>4(r)] are connected to the transverse uniaxial and biaxial director fluctuations, respectively. 
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Due to the symmetry in the (x, y) plane, it is suitable to work in a mixed (kj_, z) representation, i.e., <pi — (j>i(k±, z), 
where kx = (k x ,k y ) is the wave vector in the (x,y) plane. Now, projecting the linearized Euler-Lagrange equation 
on the base tensors g,, one obtains the following independent Schrodinger-type eigenmode equations: 



4>i(z) = Ei4>i(z), i = 0, •••,4, 



(29) 



where the functions Vo(z) = 1 + r — 6u(z) + 6u(z) 2 , V\(z) = Vz(z) = 1 + r — 3u(z) + 2u(z) 2 , and Va(z) = V±(z) 
I ■ '." Qu(z) + 2u(z) 2 may be thought of as potential energies in a related quantum-mechanical problem. The 
profile function u(z) is defined by Eqs. (|l9| ) and (po[). The correlation functions of the order-parameter fields 
Gi(k±; z, z') = (<f>i(\c±,z)<f)i(]s.j_,z')) satisfy the equations 



Gi(k±;z,z') = 8(z-z'), i = 0,---,4, 



(30) 



and can be represented, using the eigenmodes (j)[ n \z) and the eigenvalues E\' 11 of Eq. (|29|), in the form 



(n) 



(«) 



(31) 



In what follows we study the eigenvalue problem in the two extreme cases of surface free-energy potentials specified 
above. In particular, we argue that the local uniaxial director excitations are characteristic of surface potentials 
favoring the nematic state, whereas in the case of weak-anchoring potentials suppressing the nematic order, as a rule, 
the local biaxial exciations are favored. In principle, the mode spectrum may be disturbed by several issues, such as 
higher-order fluctuation effects, finite-size corrections, long-range substrate forces, and the order electricity. Some of 
the mentioned effects modifying the local director excitations are discussed below. As to the external electromagnetic 
fileds, this problem deserves a separate study. 



A. Strong-anchoring limit: local uniaxial director modes 

As discussed above, for substrates enhancing the scalar order parameter (u = u s > 1) (and for a fixed surface 
homeotropic orientation of no), the surface phase diagram contains a line of critical points describing a complete 



wetting. The related eigenvalue problem, Eq. (29), has recently been numerically analyzed by Ziherl and Zumer [ p. 2| . 
In this strong-anchoring case, the fluctuation modes are pinned at the surface, i.e., 4>i(z)\ z=0 = 0, i = 1, • • • , 4. Two 
types of local exitations, related to the potential wells shown in Fig. 2, have been established. The first one is the 
lowest soft mode 4>q (z) characterized by an energy Eq(t) cx t and related to the scalar order-parameter field: This 
mode describes fluctuations of the mean- field interface position located at z = di. Physically, this excitation appears 
as a result of the broken translational symmetry [ u(z) ^ 0] and its existence is a typical feature of the critical wetting 
transition. The second type of localized modes found in the discussed work are connected with the uniaxial director 
fluctuation fields 4>i(z) and cj>2 (z). In the vicinity of Tni the uniaxial modes soften and become gapless at Tmi- As 
seen below, these features of the local uniaxial director excitations can easily be obtained analytically in the limit 
t -► 0+ 

The physics behind the local uniaxial director modes is simple. At a mean- field level and in the limit r — + + , the 
interface width, as a rule being of the order of the correlation length £ , is much smaller than the layer thickness 
di ~ £ohi(l/r), so that the region < z < di may be thought of as a nematic plate of width di, Eq. (pl|). In the 
same limit, and for the low-energy excitations, one can use the following simpified form of the potential Vi(z) (see 
Fig. 2): Vi(z) m in the interval < z < di, and V\(z) w 1 for z > di. Thus, using the continuity of the logarithmic 

derivative of the field <f>i(z) at z = di, one finds the excitation energies E^ = £o^n- The parameter k n satisfies the 
implicit equation 

k n di = 7rn - arcsin(fc n £ ), n = 1, 2, • • • , n max , (32) 

where the number of localized modes n max is finite, fixed by the condition < fc„£o ^ 1j an d dependent on the 
reduced temperature r ( since di = di(r)). In the limit r — * + , when di ~ £ohi(l/r), the excitation energies E^ 
take the followind asymptotic form: 
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Ei — 2/i/ s ' n — 1 ^ 2, • • ■ , n m ax ■ (33) 
In (1/t) 

This expression is in agreement with the numerical results of Ref. ]l2| ] and reproduces, in particular, the observed 
cusplike bihavior of the low-laying energy levels. We see that the low-energy levels give a direct information about 
the logarithmically divergent interface position, dj ~ £ohi(l/r). This critical behavior is characteristic for the short- 
ranged substrate interactions. For example, the Van der Waals type interactions imply the power-low critical behavior 
di cx r" 1 / 3 , whereas the Freedericksz-type wetting transition was shown to give the asymptotic form dj cx t -1 / 2 [gj] . 
It is clear that the above excitation spectra can give a valuable information on both the nature of presuface forces 
and the criticality by itself. 



B. Weak-anchoring case: local biaxial director modes 

Let us remember that in the case of weak-anchoring substrate potentials the complete wetting regime is realized for 
&S > Co" 1 as t — > 0~. Now it is the isotropic phase which wets the wall and the presurfacc layer of thickness d/y — djv( r ) 
may be thought of as a plate occupied by the isotropic phase. As in the limit of strong anchoring, the critical behavior is 
controlled by the lowest scalar order-parameter mode (j>^ describing fluctuations of the mean-field interface position at 
z = dtq ■ A good approximation for (f>^ can be obtained in the limit r —> CP, when the interface profile u(z), Eq. fl2C|), 
practically coincides with the profile of the infinite system at coexistence, u(z) as u°°(z) = {1 + tanh[(z — djv)/2£o]} 
(exlcuding a small vicinity of the substrate). Since in the infinite system the function (^(z) = du°°(z)/dz is a solution 
of Eq. ( J29| ) with the potential Vq(z) = Vo[u°° (z)], it is clear that the variational ansatz 

4 0) (z) = C(r)e- Z '^ + 5 , (34) 

° ' K ' Acosh 2 [(z-d N )/2^Y V ' 

give s a good approximation for <fi^ (z) in the limit r — > Q — . Here the interface position <ijv = ^(r) is defined by Eq. 
( |22| ) . The first term in Eq. (|34| ) is introduced in order to satisfy the boundary condition ([h]) , whereas the second one 
is a solution of Eq. ( p9| ) for the infinite system at T = Tnj. From Eq. ( |To| ) one finds that the parameter C — C(t) 
takes in the limit r — > CP the asymptotic form C(r) = — (a s £ — l)w (r), where u a (r) oc \t\ x / 2 is given by Eq. fl24| ) 
in the case a s ^o > 1- Using the above ansatz as a variational function, one can show analytically p2[ that the lowest 

excited state (z) has the energy Eq ^ cx |r| and is connected to fluctuations of the mean-field interface position at 
z = dpf. More details on the importance of the discussed soft excitation will be presented in the next Subsection. 

Now, let us address the nematic-director fluctuations. Due to the fact that biaxial director fluctuations are strongly 
suppressed in the bulk nematic phase (as compared to the isotropic bulk phase), there is a well in the potential function 
Vs(z) (see Fig. 3) in the eigenmode equations (p^). The potential well is located in the presurface quasi- isotropic 
domain with a characteristic thickness d^. Thus, in the weak-anchoring case bound biaxial director modes located in 
the persurface domain can appear. On the other hand, the uniaxial director modes, being gapless Goldstone modes in 
the bulk nematic phase, are controlled by the monotonically decreasing potential V\ (z), so that they will be delocalized 
all over the sample. In the vicinity of Tmi the low-energy levels of the local biaxial modes can easily be obtained 
analytically by use of the same procedue applied above to the uniaxial excitations in the regime of strong anchoring. 
Now the potential V%{z) is simplified as V${z) m 1 (for < z < dw), and Vs(z) w 9 (for z > djv)- The boundary 
condition at z — for the biaxial field ips (z), as obtained from Eq. (|6|), reads 



#3(2) , , , 
a s (f>z{z) 



dz 



= 0. (35) 

2=0 



Using the same arguments as in the previous case, one finds the excitation energies E^ n) = 1 + ilk 2 n , where the 
parameter k n now satisfies the followig implicit equation: 

k„dM = 7m — arcsin | ° " | — arcsin | — " | , n — 1, 2, • • • , n max . (36) 

The number of localized modes n max is fixed by the condition < k n ^/2^/2 < 1, and n max also depends on the 
reduced temperature r through the function djy — dj^ir). Using the asymptotic expression djv — Col n (l/l T l), onc 
gets the following result in the limit r — >• 0~: 
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E { 3 n) = 1+ " " - , n = l,2,---,n mox . (37) 



In (1/ 

We see that in the limit r — > 0~ the local biaxial modes are stongly softened as compared to the bulk biaxial 

?3™^ change from £"3" 



fluctuations in the nematic phase: For every n, the energies E^ n) change from E ( 3 n) = 9 (the gap in the bulk nematic 



phase) to Ef> = 1 (the gap in the bulk isotropic phase). As expected, the local biaxial mode spectrum is gapped at 

C. Role of the interface-position fluctuations 

Since the upper critical spatial dimension for a critical wetting transition is D = 3 p^JT^ | , higher order fluctuation 
effects are not excluded. Indeed, as discussed above, the singular part of the mean-field free energy in the regime of 
complete wetting behaves like f* m 9 cx — rln(r), Eq. p3). On the other hand, the one-loop fluctuation contribution 

to the free energy is related to the lowest soft mode <p^\z) since the energies of the excited states of Eq. ( p9| ) are 
separated from with a finite gap. Using this fact, it is easy to see that the fluctuation part of the free energy is 
//'(t) cx —Eq ^ Ixi(Eq) ~ — rln(r), so that it compares to the singular part of the mean-field free energy. Remember 
that for uniform-director configurations the scalar order-parameter field 0o( r ) is decoupled from the director fields 
in the Gaussian approximation. Therefore, in this approximation most of the results known from the scalar theory 

p2| can be applied to nematic systems without changes. However, the fluctuation mode (fi^ will effectively disturbe 
the studied local director excitations. It is easy to see this qualitatively if remember some of the fluctuation effects 
connected to the soft mode 4>^(z), originally obtained in a context of the scalar order-parameter theory. At the first 
place, it can be shown that this mode produces the singularities in the correation function 

(to(rL,*i)0o(O,*2)) ^exp(-r ± /e||)4 0) (zi)^ 0) (^2), (38) 

where £|| = \/ \J E^ cx I/a/t is the characteristic length of the scalar order-parameter correlations parallel to the 
surface. Thus, there are critical long-range correlations of the scalar order-parameter fluctuations parallel to the 
surface (capillary waves). Another criticality connected to (ffi (z) is the predicted divergence of the interfacial width. 
Denoting by di(r±) = dj + ((r±) the local interface position, it can be shown that the characteristic interface thickness 
£± = y/ < ( 2 > diverges as r — > + according the asymptotic form £j_ oc ln£|| cx \fd~i ]l4j . Thise critical increase of 
the effective interface width, clearly, will effectively change the potential V\{z) in Eq. (|29|) which controls the director 
excitation spectrum. On the other hand, the threatment of the local fluctuation modes used throughout the paper 
remains valid, since close enough to T/vj the effective interfacial width will be much smaller than the characteristic 
domain thickness. 

D. Role of the order electricity 

There is a series of issues which might invalidate the uniform-director appoximation suggested in the above analysis 
[^J,^5|. The order electricity | f26|]2l| belongs to this series. Since the nematic- isotropic interface is characterized by a 
strong variation of the scalar order parameter, one can expect that the order-electricity effect will play an important 
role in the wetting phenomena. Let us concentrate on the strong-anchorig limit, when the nematic director is strongly 
anchored along the normal to the surface. The inhomogeneity in the profile function u(z), Eq. ([Toj) , generates along 
the axis z the electric polarization field P* — (fi^Oz + ^2) (du/dz), where n and T2 are the order-electric coefficients 
p7[ . Denoting by 9 the polar angle of the director no = (sin 9, 0, cos 9), it is easy to see that 9 = 9(z) since the 
free-energy density term f a , connected with P", mixes no and du/dz: 



1 2ir , „ . 2 / du\ 



/ = -iW = ^)(-cos^ + , 2 ) 2 ^j 2 . (39) 

Here e zz (9) = e± + (e^ — e jj cos 2 is the zz component of the dielectric tensor and E* is the z component of the 
induced electic field. Combaining the last equation with the elastic term in Eq. (||), one finds the result that the 

elastic constant L in Eq. (|J) is effectively renormalized as L — > L e ff(6) = L + [4n / e zz (6)] (ri cos 2 9 + r-i^ . To find 
the function 9{z) we will follow the arguments of Ref. p7[ |, applicable in our case for r — > + . The Euler-Lagrange 
equation for the polar angle 9, as obtained from Eq. (0), reads 



9 



dz 2 



1 dL 



eft 



6Lu 2 d9 



du 
dz 



(40) 



As shown above, du/dz oc </>o°(z) is a 5-like function centered at z = cij, Eq. (|3~4|). Thus, the function #(z) can be 
approximated with a linear function in the region < z < di, excluding the vicinity of the interface and, possibly, of 

1/2 

the surface. Since the interfacial energy is proportional to L e ff , the minimum condition for this energy dL e ff/d9 = 



also defines the average polar angle at the interface 9q 28 1 . We do not consider here possible deviations from the linear 
behavior of 6(z) close to the substrate 
by the polar angle 



29[. Therefore, as t — ► + the inhomogeneous director state is characterized 



9(z) = 9 —~9, 



1 



T 



£o ln(l/r)' 



H 



z < dj. 



(41) 



The inhomogeneous director state described above can generate several changes in the fluctuation mode dynamics. 
Since the local base tensors are now coordinate-dependent, gi — gi(z), the normal mode equations ( |29] ) are coupled. 
In particular, there is a coupling between the longitudinal scalar order-parameter field <fio(z) and the transverse 
director fluctuations 4>2{z) and 4>s{z). The mode coupling is, however, asymptotically small, as it is controlled by the 
small parameter £o/di oc ln -1 (l/r). As the interaction decreases logarithmically , the effect of order-electicity mode 
coupling can play an important role in real experiments. 



IV. SUMMARY OF THE RESULTS AND DISCUSSION 



In this article the emphasis was on specific features of the presurface fluctuation mode dynamics in nematic liquid 
crystals, arising in a regime of critical wetting when the system is close to the nematic-isotropic phase transition 
temperature. The importance of substratate anchoring for the director mode dynamics is demonstrated by use of two 
extreme types of surface free-energy potentials. 

Our key results are as follows: In the case of uniform-director mean- field configurations and for a homeotropic 
director anchoring, surface free-energy potentials strongly suppressing the orientational order favor the appearence 
of bound biaxial nematic-director fluctuation modes located in the domain occupied by the wetting thermodynamic 
phase. Instead, substrate forces which enhance the orientational order promote bound uniaxial modes in the presurface 
domain. Close to the phase coexistence temperature both types of local director excitations are strongly softened 
(as compared to their bulk counterparts) and acquire characteristic cusplike low-energy spectra which are directly 
connected to the critical behavior of the mean-field interface position. Since the physics of these local excitations 
closely reflects the wetting critical behavior, they can be used as a tool to examine the critical wetting phenomena 
and the presurface interactions in nematic liquid crystals. 

We have considered one of the possible reasons for the experimentally observed different spectra of elementary 
exitations, ranging from very small to very large (as compared to the bulk) relaxation rates. Clearly, there could be a 
number of other factors modifying substantially the director exitation spectra. We have already briefly discussed two 
of them, i.e., the higher-order fluctuation effects and order electricity effects. As to the interface position fluctuations, 
in the upper spatial critical dimension D = 3 they lead to a divergence of the effective interface width. Thus, the 
local excitation specta can be disturbed, as the potentials Vi(z) will be effectively changed as well. In addition, the 
fluctuation effects in D = 3 are known |Q to depend on the value of the dimensionless parameter uj — T/An^a, 
where £ is the bulk correlation length of the phase attracted to the wall and a is the surface tension of the free 
interface. It will be interesting to see the role of this parameter in different liquid crystal materials. As to the order- 
electricity effects, they are expected to play an important role in the real experiments due to the strong variation of the 
scalar order parameter in the interfacial region. The coupling of the normal Gaussian modes is the main effect of the 
inhomogeneity. It is important that in this case the director fluctuation modes are coupled with the " dangerous" scalar 
order-parameter soft mode, so that the picture of the Gaussian fluctuation dynamics can be importantly disturbed. 

In conclusion, the finite-size effects, the long-ranged substrate interactions, and the external electromagnetic fields 
belong to the class of factors that will modify the discussed local fluctuation mode dynamics. Clearly, these issues 
deserve future studies. 
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FIG. 1. Scalar order-parameter profiles for the reduced temperatures r = —10 2 , —10 4 , and —10 6 in the case of substrates 
strongly suppressing the nematic order. The inflaction point at z = cLm marks the characteristic isotropic domain thickness. 
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FIG. 2. The fluctuation mode potentials in Eq. (ESJ) in the srtong-anchoring limit. The wells in Vo{z) and Vi(z) are connected 
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